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Abstract. This is a survey on discrete linear operators which, besides approximating in 
Jackson or near-best order, possess some interpolatory property at some nodes. Such operators 
can be useful in numerical analysis. 



A central problem in approximation theory is the construction of simple functions that ap- 
proximate well a given set of functions. Traditionally, by "simple" we mean polynomials or rational 
functions as they are easily implemented on computers, and the set of functions is generally char- 
acterized by continuity or belonging to some LP class. These functions may be defined on finite or 
infinite intervals, or in complex domains. Another characteristic is the measure of error which can 
be the supremum or LP norm, etc., with a weight usually introduced in the case of non-compact 
domains. 

But what are the available data to construct such approximations? From a practical point 
of view, this is a crucial question. For example, convolution integrals are very useful in proving 
Jackson's theorem, but to actually calculate them one needs complete knowledge of the function. 
Obviously, this precludes any numerical application. In practice, we are generally given a discrete 
set of data (like function values at certain points) from which we wish to reconstruct the function. 

Another feature which is important to us is that we prefer to construct linear operators since 
they are easier to handle. For example, the so-called "best approximation" is not at all "best" for 
our purposes, since (except in inner product spaces) it represents a non-linear operator which may 
be very difficult to calculate. 

A natural candidate that satisfies the above requirements is some kind of interpolation operator. 
It is based on a discrete set of data, and additionally gives a zero error at an increasing number of 
well defined points. Depending on the nature of the interpolation, it is relatively easy to construct 
such interpolation operators. 

However, there is a significant drawback to these operators. In the case of Lagrange interpo- 
lation, for example, no matter how we choose the n nodes we get at least an extra O(logn) factor 
compared to the optimal Jackson order of convergence. And if we try to avoid this problematic 
situation by considering Hermite-Fejer interpolation then, although we obtain the Jackson order 
of convergence in some cases, this process is saturated. That is, it will not give the Jackson order 
of approximation, for example, for Lip 1 classes, let alone for classes of functions with higher order 
of smoothness. 
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It transpires that the reason for this negative behavior of interpolating polynomials is the strict 
restriction on the degree. In 1963, at an Oberwolfach conference, Paul Butzer raised the following 
question: Is it possible to construct an interpolation process that gives the Jackson order of conver- 
gence? Of course, this was meant in the sense that we should now allow interpolating polynomials 
of degree higher than Lagrange or Hermite-Fejer interpolation, say polynomials of degree < an for 
some constant a. It was Geza Freud [11] who answered this question affirmatively by constructing 
polynomials p n of degree at most 4n — 3, interpolating at ~ n/3 nodes and approximating to the 
Jackson order 

max \f(x) - p n {x)\ < cu(f, 1/n), \x\ < 1, 

|x|<l 

where u is the ordinary modulus of continuity This was later improved upon by Freud and Vertesi 
[14] who constructed a sequence of discrete linear polynomial operators (DLPO, in what follows) 
J n (f,x) of degree at most 4n — 2 which interpolate at the Chebyshev nodes cos((2fc — l)/(2n))7r, 
k = 1, . . . , n, and provide a Timan type (pointwise) estimate 

- J n (f, x)\ < c[u(f, ^-^-) + i)], \x\ < 1. (1) 

Freud and Sharma [12, 13] further extended this result to operators based on general Jacobi nodes, 
and also succeeded in decreasing the degree of the polynomial to n(l + e), for an arbitrary e > 0. 
(Of course, the c in the above then depends upon e.) 

Freud's work [11] initiated a substantial series of papers exhibiting constructions of a similar 
character. R. B. Saxena [19] succeeded in constructing a sequence of DLPOs J* which realized the 
even stronger Telyakovskii-Gopengauz estimate 



\f(x) - r n (f,x)\ < cu(f, ^-^), \x\ < 1. 

Furthermore, for 27r-periodic continuous functions, O. Kis and P. Vertesi [17] constructed a very 
simple interpolating process: Let t^{x) be the fundamental functions of trigonometric interpolation 
based on the equidistant nodes x^ = 2kir/ (2n+l), k = 0, . . . , 2n, i.e., let £k{%) be that trigonometric 
polynomial of order n for which £k{xj) = Sjk, j, k = 0, . . . , 2n, and consider the sequence of operators 

n 

U n (f,x) = ^/(x fc )[44(x) 3 - 34(x) 4 ]. 

fc=0 

Then, for every 27r-periodic continuous functions /, we have 

max|/(x) - U n (f,x)\ < cw(/, -), 

x€R n 

and, for this trigonometric polynomial U n (f,x) of degree at most An, evidently 

U n {f,x k ) = f(x k ), k = 0, ...,2n, 

because of the nice identity 

n 

^[44(x) 3 -34(x) 4 ] = l 

k=0 
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(cf. Turetskii [27]). With a proper transformation, Kis and Vertesi [17] also constructed an operator 
that realized the Timan estimate (cf. (1)) and interpolates at some nodes. In that construction, 
the limit of the ratio of the degree of the polynomial and the number of points of interpolation 
was greater than 1. Later, O. Kis and J. Szabados [16] sharpened this to the following: Given 
< e < 1, let n > 20/ e 2 and f G C[— 1, 1]. There exists a sequence of DLPOs p n of degree at most 
n(l + e) such that p n interpolates f in at least n points and 



\f(x) - Pn (x)\ < ^ U(f, ^ ), |*| < 1. 

This answers a question raised by Freud and Sharma [12] about the construction of linear operators 
of minimal degree compared to the number of points of interpolation, that at the same time realize 
the Telyakovskii-Gopengauz estimate. 

In the algebraic case, the above mentioned operators are based on nodes of orthogonal polyno- 
mials which is not the best choice from a numerical point of view. Usually, data about an unknown 
function are given at equidistant nodes. Since Lagrange interpolation at equidistant nodes behaves 
very poorly (the norm of the operator grows exponentially with the number of nodes), it was not at 
all obvious if convergent operators based on equidistant nodes could be constructed. On p. 581 of 
[3], R. DeVore raised this question while requesting a Jackson order of convergence. One possible 
solution is the following: As an intermediate approximation, we take the rational functions 

^— ' (x — k/nr / ^— ' Ix — k/nr 

k=—n I k=—n 

based on equidistant nodes. These are a special case of the so-called Shepard operators which 
approximate to Jackson order (for details, see later). Then, to approximate (2) by polynomials, 
we can use any standard constructive operator (for example, the one developed by Bojanic and 
DeVore [2]). An easy calculation shows that the moduli of continuity of / and (2) are equivalent, 
thus solving the problem (for details, see [24]). Also, if the operator obtained in this way is of 
degree n, then for some c, < c < 1, it is possible to modify it so that the new operator will 
interpolate at the equidistant nodes k/(cn), k = 0, ±1, . . . , ±cn (cf. [25], Theorem 3). 

With the exception of Lagrange interpolation, all of the above mentioned operators are satu- 
rated, i.e., the order of convergence cannot be improved upon beyond a certain limit, even if the 
function has better and better structural properties. The question then becomes: Can we con- 
struct discrete linear operators that are not saturated, and approximate "close" to the order of 
best approximation E n (f)7 Of course, by Korovkin's theorem, such operators cannot be positive. 
In the case of 2-7r-periodic functions, the well-known de la Vallee Poussin means do have the de- 
sired approximation order, but they are not discrete. In [23], the following family of operators was 
constructed. Let j, n be positive and k, I nonnegative integers such that \{jn + kn — k + £ — 1) is 
a nonnegative integer. Set 

_ sin^(sinf ) fc (cosf)^ 
S3kUt) ~ j(nsin|)*+i ' 

2-kv 

t v = —, i/ = 0,±l,±2,..., 
jn 

and let 
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jn—l 

Sjk£n(g,t) = 9(tv)Sjken(t), 
u=0 

be a discrete linear operator denned on all 27r-periodic continuous functions g. This is a gener- 
alization of many of the classical kernels and operators, sioon(^) ( n odd) is the Dirichlet kernel, 
snon(t) is the Fejer kernel, ssio n (t) is the de la Vallee Poussin kernel, and si30n(i) is the Jackson 
kernel. Correspondingly, Sioo n (g,t) [n odd) and S\o\ n {g,t) (n even) are ordinary interpolation 
polynomials, Sno n (g,t) is the Jackson polynomial and S3io n (g, t) is the discrete version of the de 
la Vallee Poussin operator (see [16] and [23]). These sequences of operators converge if the so-called 
Lebesgue constant 

jn—l 

\\Ljkin\\ = max ^ \Sjkin{t - t v )\ 

e I/=0 

remains bounded as n — > oo (the other parameters being considered as constants). Namely, the 
following typical error estimate holds: 

m.ax\S jMn (g,t) - g(t)\ < (1 + \\L jken \\)E^(g), (3) 

t GR. 

where £'„(•) is the error in the best approximation by trigonometric polynomials. For some special 
choices of the parameters j, k, £, the following relations hold: 



3 



1 ^ 2u - 1 2„ 8?r 

\\LjWn\\ = ~ Z^ cot A ■ 11 = -( lo g— +7) + « 
J ~f 4j IT J 

where 7 = 0.5772 • • • is the Euler constant and < ctj < 

2 

||-^jiin|| < — log j + 2.283 if j is even, n is odd; 

7T 

\\L22mW < t|; 

1 1^33271 1 1 < 

In all of these cases, (3) yields the order of best approximation. In addition, evidently 

Sjk£n(g,tu) = g(t„), v = 1,... ,jn, 
i.e., the ratio of the order of the operator to the number of interpolation points is 

jn + kn + £ — k + 1 
2jn 

which is fairly close to the optimal ratio 1/2 (which corresponds to the minimal degree interpolation) 
provided j is large compared with k. 

The obvious transformation x = cost converts these results into estimates for continuous 
functions on the interval [—1, 1]. Let 
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Pjk£n{f,x) = S j k e n (f (COS t),t) 

and 

x v = cost v , v = 0, . . . , jn — 1. 
This new polynomial operator will have the analogous properties 

Pjkin(f,x v ) = f(x v ), v = 0, jn- 1, 

and 

max|/(ar) - P jktn (f,x)\ < (1 + \\L jktn \\)E n (f), 

\x\<l 

where E n (-) is the error in the best approximation by algebraic polynomials of degree at most n. 

Let the number r n = m n /n be the ratio of the degree m n of a polynomial operator to the 
number n of nodes where it interpolates. For Lagrange interpolation we always have lim^^oo r n = 1. 
In the case of (convergent) Hermite-Fejer interpolation, lim^oo r n = 2. A classical result of 
Bernstein [1] says that if e > is arbitrary, then there exists a linear polynomial operator with 
bounded norm such that 



lim r n = 1 + e. (4) 

n — >oo 

In general, it is desirable to investigate the relation between the norm ||L re || of a (not necessar- 
ily linear) polynomial operator of degree m n and the number of nodes n where it interpolates. 
Generalizing the results just mentioned, it can be shown that for such operators 

ll^nll 

lim - — ^ > 0. 

ra^oo log — -7; 

This was proved first for Chebyshev nodes [26], and later generalized by B. Shekhtman [20] to an 
arbitrary system of nodes. We emphasize that this result is true without assuming linearity of the 
operator. It holds, for example, for operators of the form 

n 

L n (f,x) =Pn(f,x) + ^2\f(x k ) -Pn(f,X k )]q k (x), 
k=l 

where p n is the best approximating polynomial of degree m n and the q k are polynomials of the 
same degree such that q k (xj) = S k j, j, k = 1, . . . , n. Sequences of operators having property (4) for 
general systems of nodes x k = cos t k , k = 1, . . . , n, and approximating to the order 0(-E n (i+ £ )) were 
constructed by Erdos, Kroo and Szabados [10] who showed that such approximation is possible if 
and only if 

r NJJn).l , r in 

hmsup — < — whenever iim n\l n \ = oo, 

n^oo n\I n \ IT n^oo 

and 



lim n ■ min (t k — t k+ i) > 0, 

n^oo l<k<n— 1 
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where N n (I n ) is the number of tk's in the interval I n C [0, ir]. 

We have, until now, dealt with polynomial operators. Prom a numerical point of view, rational 
functions are just as useful as polynomials, and as we shall see, they enjoy certain advantages over 
polynomial operators. We have already briefly mentioned the operator (2). With an exponent 
2 instead of 4, and in the case of bivariate functions, this operator was originally introduced 
by D. Shepard [21]. His interesting paper went unnoticed for some time, but his operator was 
rediscovered in 1976 by J. Balazs (oral communication). It was then proved in [24] that for the 
sequence of operators 

R n (f,x) = Y , /(fc /" ) , / Y - t ]~r-r , (5) 

w ' 1 ^ x-k h\p / ^ x-k nP v ' 

k=—n I k=—n 

in the case p = 4, we have 

max \f(x) - R n (f,x)\ < cu(f, -), 
|x|<i n 

i.e., a Jackson order of convergence. The major advantage of this sequence of discrete linear opera- 
tors is that it is based on equidistant nodes (and not on the roots of some orthogonal polynomials). 
In addition, these operators (for p = 2) have the properties 



R n (f,k/n) = f(k/n), R' n (f,k/n) = 0, k = 0,...,±n, 

i.e., they behave like Hermite-Fejer interpolating polynomials. 

Balazs' rediscovery of Shepard's operators initiated a series of papers dealing with further 
problems and generalizations. We mention only a few significant papers and results. G. Somorjai 
[22] proved that, for p > 2, 

max|/(x) -R n (f,x)\ < - & /eLipl, 
| cc | < l n 

and 

max \f(x) - R n (f,x)\ = o (J^j 4^ f = const., 

i.e., the saturation property of the operator R n . The important (and much more difficult) case 
p = 2 was handled by Delia Vecchia, Mastroianni and Totik [4]. They proved that for p = 2 the 
trivial class (for which the approximation order is o(l/n)) is the set of constant functions, and the 
saturation class (for which the order of approximation is 0(l/n)) is contained in n Q <iLipa. They 
also proved the surprising fact that, for arbitrary subsequences of the Shepard operators, arbitrarily 
fast convergence may be achieved for some nonconstant functions. 

The case < p < 2 was handled by X. L. Zhou [28]. With the notation 

T e, P f ■= ma 

0<x 

he proved that the saturation order is n 1_p or 1/logn, and the saturation class is 



I 



/w - / w dt 



0<t<l, |t-x|>£ \t — x\ p 



{/ : / G Lip (p - 1), sup T £jP f < cx)} 

£>0 
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or 

{/ : u(f,t) = 0{\ \ogt\~ 1 ), sup T eA f < 00} 

£>0 

according asl<p<2orp = l, respectively. He also showed that the case < p < 1 is not 
interesting, since then 

max \f(x) — R n (f, x)\ — > 44> / = const. 

0<a;<l 

Many of the above mentioned results have been generalized to weighted approximation (see 
[6]-[9], [18]), as well as to general systems of nodes (see [5]). 
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